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Short wave equations were introduced in connection with the nonhnear reflection of weak 
shock waves. They also relate to the modulation of a gas-fluid mixture. Khokhlov-Zabolotskaya 
i equation are used to describe the propagation of a diffraction sound beam in a nonlinear medium. 

We give a new algebraic method of solving these equations by using certain finite-dimensional 
stable range of the nonhnear terms and obtain large families of new explicit exact solutions 
J> ' parameterized by several functions for them. These parameter functions enable one to find the 

(SJ ! solutions of some related practical models and boundary value problems. 

cn 

^ ; 1 Introduction 

> : 

■ Khristianovich and Rizhov [8] (1958) discovered the equations of short waves in connection 
?H ' with the nonlinear reflection of weak shock waves. The equations are mathematically 
equivalent to the following equation of their potential function u for the velocity vector: 

2Utx - 2(X + + Uyy + 2kur, = 0, (1.1) 

where k is a real constant. For convenience, we call the above equation "the short wave 
equation". The symmetry group and conservation laws of (1.1) were first studied by 
Kucharczyk [14] (1965) and later by Khamitova [6] (1982). Bagdoev and Petrosyan [2] 
(1985) showed that the modulation equation of a gas-fluid mixture coincides in main 
orders with the corresponding short-wave equation. Roy, Roy and De [23] (1988) found 
a loop algebra in the Lie symmetries for the short-wave equation. Kraenkel, Manna and 
Merle [13] (2000) studied nonlinear short-wave propagation in ferrites and Ermakov [3] 
(2006) investigated short-wave interaction in film slicks. 
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Khokhlov and Zabolotskaya [7] (1969) found the equation 



2Utx + iuUx)x - Uyy = 0. 



(1.2) 



for quasi-plane waves in nonlinear acoustics of bounded bundles. More specifically, the 
equation describes the propagation of a diffraction sound beam in a nonlinear medium 
(cf. [4], [20]). Kupershmidt [15] (1994) constructed a geometric Hamiltonian form for the 
Khokhlov-Zabolotskaya equation (1.2). Certain group- invariant solutions of (1.2) were 
found by Korsunskii [12] (1991), and by Lin and Zhang [16] (1995). The three-dimensional 
generalization 



and its symmetries were studied by Krasil'shchik, Lychagin and Vinogradov [17] (1986) 
and by Schwarz [25] (1987). Martinez-Moras and Ramos [18] (1993) showed that the 
higher dimensional classical W-algebras arc the Poisson structures associated with a 
higher dimensional version of the Khokhlov-Zabolotskaya hierarchy. Kacdryavtsev and 
Sapozknikov [9] (1998) found the symmetries for a generalized Khokhlov-Zabolotskaya 
equation. Sanchez [24] (2007) studied long waves in ferromagnetic media via Khokhlov- 
Zabolotskaya equation. Morozov [19] (2008) derived two non-equivalent coverings for 
the modified Khokhlov-Zabolotskaya equation from Maurer-Cartan forms of its sym- 
metry pseudo-group. Rozanova [21, 22] (2007, 2008) studied closely related Khokhlov- 
Zabolotskaya-Kuzentsov equation from analytic point of view. Kostin and Panasenko [11] 
(2008) investigated nonlinear acoustics in heterogeneous media via Khokhlov- Zabolotskaya- 
Kuzentsov-type equation. All the above equations are similar nonlinear algebraic partial 
differential equations. 

Observe that the nonlinear terms in the above equations keep some finite-dimensional 
polynomial space in x stable. In this paper, we present a new algebraic method of solving 
these equations by using this stability. We obtain a family of solutions of the equation 

(1.1) with k = 1/2,2, which blow up on a moving line y = f{t). They may refiect 
partial phenomena of gust. Moreover, we obtain another family of smooth solutions 
parameterized by six smooth functions of t for any k. Similar results for the equation 

(1.2) are also given. Furthermore, we find a family of solutions of the equation (1.3) 
blowing up on a rotating and translating plane cosa(t) y + sinQ;(t) z = f{t), which may 
reflect partial phenomena of sound shock, and a family of solutions parameterized by 
time-dependent harmonic functions in y and z, whose special cases are smooth solutions. 
Since our solutions contain parameter functions, they can be used to solve certain related 
practical models and boundary- value problems for these equations. 

On the list of the Lie point symmetries of the equation (1.1) in the works of Kucharczyk 
[14] and of Khamitova [6] (e.g. cf. Page 301 in [5]), the most sophisticated ones are those 
with respect to the following vector flelds: 





a'ydx + ady+ xy{a" + a') - ^{a'" + {k + l)a" + ka') 9„, (1.4) 



X2 = Pd, + [y\p" +{k + 1)/?' + kp) - x{(5' + /?)]a„, 



(1.5) 



where a and f3 are arbitrary functions of t. Among the known Lie point symmetries of 
the Khokhlov-Zabolotskaya equation (1.2) in the works of Vinogradov and Vorob'ev [26], 
and of Schwarz [25] (e.g. cf. Page 299 in [5]), the most interesting ones are those with 
respect to the following vector fields: 

1 1 

^3 = ^oi'yd:^ + ady - -a"ydu, (1.6) 

X4 = f3dt + + -13 ydy du- (1.7) 

b 6 6 

The symmetries of the three-dimensional Khokhlov-Zabolotskaya equation (1.3) causing 
our attention are those with respect to the vector fields (e.g. cf. Page 301 in [5]): 

= IQt^dt + (4te + 3y^ + 3z^)9^ + mydy + mzd;, - {Ax + mu)du, (1.7) 

1 1 

= i^oc'yd^ + adz - ^a"ydu, (1-8) 

X7 = ^P'zd, + Pdy - ^I3"zd^. (1.9) 

We find that the group-invariant solutions with respect to the above vector fields Xi- 
X'j are polynomial in x. This motivates us to find more exact solutions of the equations 
(1.1)-(1.3) polynomial in x. 

In Section 2, we solve the short-wave equation (1.1). Although the equation (1.2) 
can be viewed as a special case of the equation (1.3), we first solve (1.2) in Section 3 for 
simplicity because our approach to (1.3) involves time-dependent harmonic functions and 
sophisticated integrals. The exact solutions of the equation (1.3) will be given Section 4. 

2 Short Wave Equation 

In this section, we study solutions polynomial in x for the short wave equation (1.1). By 
comparing the terms of highest degree in x, we find that such a solution must be of the 
form: 

n = f{t, y) + g{t, y)x + h{t, y)x^ + C{t, y)x\ (2.1) 

where f{t,y), g{t,y), h{t,y) and ^(t,y) are suitably-differentiable functions to be deter- 
mined. Note 

u^^g + 2hx + S^x'^, u^^ ^2h + 6^x, (2.2) 

UtX = 9t + '2htX + 3^tX^, Uyy = fyy + QyyX + kyyX^ + ^yyX^ , (2.3) 

Now (1.1) becomes 

2{gt + 2htx + 3^tx^) - 2{g + {2h + l)x + 3^a;2)(2/i + 6^a;) 

+fm + 9yyX + hyyX^ + iyyX^ + 2 A) ( ^ + 2hx + ?>ix^) = 0, (2.4) 



which is equivalent to the following systems of partial differential equations: 

^yy — 36^ , 



(2.5) 



hyy = 6^{6h + 2-k)- (2-6) 

Qyy = + 4(1 - k)h + I2ig - 4ht, (2.7) 

fyy = Agh - 2gt - 2kg. (2.8) 

First we observe that 

C = r- ^ (2.9) 

{V6y + P{tW 

is a solution of the equation (2.5) for any differentiable function P of t. Substituting (2.9) 
into (2.6), we get 

Denote by Z the ring of integers. Write 

h{t,y) = J2<t)i^y + Pit)y- (2-11) 

Then 

hyy = 5^ 6(z + 2){t + l)a,+2(t)(\/6y + /3(t))\ (2.12) 
Substituting (2.11) and (2.12) into (2.10), we obtain 

g + 2)(, + 1) - 6la.«(0(V6. + PWr - . ^^f^. (2.13) 

So 

-24a_i(i) = 12/3'(t), -36ao{t) = 6{2 - k) (2.14) 

and 

6(i + 4)(i-l)ai+2(t) =0, t^-2,-3. (2.15) 

Thus 

where a and 7 are arbitrary differentiable functions of t. 
Note 

and 

^2 _ 3{(3'y + A{k - 2)a (2 - fc)/j^ (fc - 2 f 

"(y6i/ + /3)4 (^|/ + /3)3 12(76?/ + /?)2 6(76?/ + /3) 36 

+2a7(V6y + /?) - /3'7(\/6y + + ^tz^l^^y + (^f + ^2^^^ + /3)6. (2.18) 

Substituting the above two equations into (2.7), we have: 

12^ A[{k + l)a + 3a'] 2{{k + 1)13' + 3p") 



9yy 



{VQy + PY {^/Qy + PY 3{VQy + l3Y 3(V6y + /?) 



Write 



- ^^^ + 16a7(V6y + /?) - 20/?'7( Vet/ + 

.!l(^±ill±M(V6, + «= + 8T^(V6, + ffl«. (2.19) 

s(t,»)=^(>i(i)(^» + /3)'. (2.20) 



Then 

^ 6[(z + 2)(z + 1) - 2]6,+2(i)(A/6?/ + (3)' 

8a^ 4p; + l)a + 3a'] 2[{k + + Sp"] 

-'-^^^±^^±^iV6y + Pf + 8f (V6, + Pf. (2.21) 

Comparing the constant terms, we get k = 1/2, 2. Moreover, the coefficients of the other 
terms give 

a' {k + l)a + 3a' {k + 1)P' + Sp" 
&-2-y, feo- g , 6i- ^ , (2.22) 

, _ 2a7 , _ P'l , _ (A; + 1)7 + _ 27^ 

&3-^, ^4-- — , 65- , bs-— (2.23) 

and 

(i + 3)i6i+2 = for i -4, -2, -1, 1, 2, 3, 6. (2.24) 

Therefore 

>,2 
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^ g ^ (fc + l)a + 3a- + 1)/^^ + 3/3^^ ^/^^ , 



3{VQy + py VQy + P 9 18 

+p(V6y + py + ^{V6y + Pf - ^(V6y + Pf 

_ (/. + l)7 + 3y ^^^ + /?)^ + |^( + /?)«, (2.25) 

where a and p are arbitrary differentiable functions of t. 
Observe that 

2a'^P' {2aa'-3aP') a' {k + l){2a' - {p')"^) 

9t — — — , + /T- — , H — — : + 



3(V6y + /3)3 3(^/6?/ + /3)2 V^y + P 18 



2q; 
H 



6 18 

2q;V 

4 



2„y + 2„-T-4(/^-)^^^^ + _ 4W±5(^+2MW^(^^ ^ 

3 81 

(Jl±M±^(^^ + ^). + ^(V6. + ^'q^i^^y^ Pf, (2.26) 



81 V - ' / 81 ' ^ 81 



^ " 3{VQy + P)' ^ Q{VQy + ^'f ^ 6{VQy + 

3{k - 2)(7 - 3aP" - 2{k + l)ap' - 3a' p' {k + + 3p'p" 

^ 18(v^y + /3) ^""^^ 36 

+ ^4 + [« 7 - ^ ] {V6y + P) 



a 

4 



{k-2)p + 6^a-4aP'j y- , {PTi {17k - 10)aj - 3a^' 

6 y P) [ Q g-^ 



3 jvv « . r-y . 162 

6 



28/3'7% (A; -2)72 3 



^(v^y + + ^^^(^1/ + + |r(v^l/ + (2-27) 
Substituting (2.25), (2.26) and (2.27) into (2.8), we obtain 

~ 3(7iy+^ 9(V6z/ + /3) 

12a(T + 2a^p' 2[p' p" - {k + l)a' - a"] 2{k + l){p'f 

^ 3iVQy + pr ^ ""^^ 3 + 9 



27 ■ L / / r ■ 3 27 



Thus 



/ = ^ + ^ — + 9 + '&y + 2apy^ + ^ lltL^y^ 

^ 27(v^i/ + /3)2 is{^y + p) y f^y 9 ^ 

Jik + l)a + 3aP" + 2ik + l)aP' + 3a'P' + 9a'^^^^ ^ ^^^^^^^^ ^ _ 

27 

^ (A: + l)a^-a^' ^2 _ 2{k + lYa ^^ ^ _ ^ + + + 1)/^" 



3 ^ 27 " ^ 9 6 108 
+^^] iVey + pr + - i^±^l|±^]( V6. + pr 

- ^^^^^^±^](V6. + ,)^ + + 1)^-7 + 3/3Y (^^ ^ 

TP (^+l£7 ^ + + ^ + ^(V6, + /3)« 

^63 15309 3402 243 ^ ^ '^^ 

-^iV6y + Pf - 27((^ + l)7 + 3y)(^ ^10 + ^^^(^6^ + /3)^^ (2.29) 
where ^ and 'd are arbitrary functions of t. 



Theorem 2.1. When k = 1/2, 2, we have the following solution of the equation (1.1) 
blowing up on the surface \/6y + i3{t) = 0: 

^3{^/6y + (3)^ V6y + P 9 18 wy^A^; 

+p{V6y + Pr +^-^iV6y + - ^{V6y + - ^^^^^^{^{V6y + Pf 
272^ 6aa + a'^p' , ik + l)iP'? 2 

_6(k + 1). + 3./>-' + 2(k + l)a/?- + 3a-/>- + 9.'^^^ ^ ^^^^^^^ ^ _ 

27 

- (A; + 1)q;' - a" 2 2(A; + l)2a ^ a^^ p'p f3'" + {k + l)p" 



+— y - y + [^- ^ + 



27 " ^ 9 6 108 

where a, P,^,a, p,9 and •& are arbitrary functions oft, whose derivatives appeared in the 
above exist in a certain open set ofR. 

When a = j = a = p = 9 = i^ = 0, the above solution becomes 

^ * - 2 P' , 2 (k+l)P' + 3P" , ^ 

"=(71^ + 1— -5(7^"" rs — (^f+'^>^ 

+ ^Ji±mly^ + . 1 ^"' " + ^1(V6. + pr. (2.31) 

Take the trivial solution ^ = of (2.5), which is the only solution polynomial in y. 
Then (2.6) and (2.7) become 

hyy = 0, gyy = 8/^2 + 4(1 - k)h - Aht, (2.32) 

Thus 

h = a{t) + p{t)y. (2.33) 

Hence 

Qyy = 4(2^2 + (1 _ k)a - a') + A{Aap + (1 - k)p - p')y + 8/3V- (2-34) 

So 

g = ^ + ay + 2{2o? + (1 - A;)a - ol)y' + ^(4a/3 + (1 - A;)/3 - P')y'' + ^/^^y^ (2.35) 



where 7 and a are arbitrary functions of t. Now (2.8) yields 

= 4(Q; + /3?/)[7 + (7?/ + 2(2a2 + (l-A;)Q;-Q;')?/' + ^(4a/3 + (l-A;)/3-/3')y' 
+ - 2[7' + o'y + 2(4aa' + (1 - k)ol - a")y^ + \{Aa' 13 + Aap' 

+(1 - - + - 2A;[7 + ay + 2{2a^ + (1 - - a')y^ 

+^(4a/? +(1- A;)/? + ^/?V] 
= 4q;7 - 27' - 2A;7 + 2{2a(T + 2/37 - cr' - k-i)y + 8((1 - A;)q; - a')(5y^ 
+4(4q;^ + 2(1 - 2k)a^ - 6aa' + k{k - l)a + {2k - l)a' + a" + (3a)y'^ 
+\{2Qa^l3 + 2(1 - ?,k)al3 - 6a(3' - Aa'(3 + {2k - l)/5' + /?" - k{l - k)i3)y^ 

+^{10aP' + (2 - 3A;)/3^ - 4/3/3')y' + ^/SV- (2.36) 

Therefore, 

, , , 9 2q;(7 + 2/?7 - cj' - A;7 o 2((1 - A;)q; - a') . 
/ = (2a7 - 7' - A;7)2/2 + ^11 ^2/' + ^^ f + r + py 



+-(4q;^ + 2(1 - 2k)a'^ - 6aa' + k{k - l)a + {2k - l)a' + a" + (3a)y^ 
3 

+^{20a^p + 2(1 - 3k)aP - Qa(3' - Aa' (3 + {2k - l)(3' + (3" - k{l - k)(3)y^ 
15 

2 AB^ 
+ -{lQa(3' + (2 - Zk)(3' - A(3(3')y' + -^y\ (2.37) 



Theorem 2.2. The following is a solution of the equation (1-1): 

u = {a + (3y)x^ + [7 + ay + 2{2a^ + (1 - k)a - a')y^ + \{Aai3 + (1 - k){i - p')y^ 

2^2 41 / 7 N 2 2q;(7 + 2/?7 - (j' - A;7 3 2((1 - A;)q; - a') , 

6 3 5 

+r + py + ^Aa^ + 2(1 - 2k)a^ - 6aa' + k{k - l)a + {2k - l)a' + a" + pa)y^ 

+^{20a^P + 2(1 - 3k)aP - Qa(3' - Aa (3 + {2k - l)(3' + (3" - k{l - k)(3)y^ 

2 AB^ 
+-i(10a/3^ + (2 - 3k)(3' - A(3(3')y' + (2.38) 
45 53 

where a, f3,j,a, p and t are arbitrary functions of t, whose derivatives appeared in the 
above exist in a certain open set ofR. Moreover, any solution polynomial in x and y of 
(1.1) must be of the above form. The above solution is smooth (analytic) if all a, (3, j, a, p 
and T are smooth (analytic) functions oft. 

Remcirk 2.3. In addition to the nonzero solution (2.9) of the equation (2.5), the other 
nonzero solutions are of the form 



e = p.(V^y + /3(t)), (2.39) 



where Pc{w) is the Weierstrass's elhptic function such that 



(2.40) 



and i is a nonzero constant and f3 is any function of t. When f3 is not a constant, the 
solutions of (2.6)- (2.8) are extremely complicated. If f3 is constant, we can take /5 = by 
adjusting l. Any solution of (2.6)-(2.8) with h ^ is also very complicated. Thus the 
only simple solution of the equation (1.1) in this case is 



u 



p,{V6y) x\ (2.41) 



3 2-D Khokhlov-Zabolotskaya Equation 

The solution of the equation (1.2) polynomial in x must be of the form 

n = /(t, y) + g{t, y)x + i{t, y)x\ (3.1) 

Then 

U^ = g + 2^X, Utx = 2^tX, Uyy = fyy + gyyX + ^J/J/X^ (3.2) 

(uu,), = d,{fg + {g' + 2fi)x + 3^^^^ ^ 2^2^3) ^ ^2 ^ ^ g^^^ ^ g^2^2_ (3 3) 
Substituting them into (1.2), we get 

2{gt + 2itx) +g^ + 2fi + Qgix + 6f - - gyyX - iyyX^ = 0, (3.4) 

equivalently, 

iyy = 6f , (3.5) 

9yy - Qg^ = 46, (3.6) 

f^^-2f^^2gt + g'. (3.7) 

First we observe that 

is a solution of the equation (3.5) for any differentiable function /3 of t. Substituting (3.8) 
into (3.6), we obtain 

Write 

^(t,2/) = 5^a,(t)(2/ + /5(t))\ (3.10) 



Then (3.9) becomes 

+ 2)(z + 1) - Q]aMt)iy + mY = 'J^f^s' (^.11) 



Thus 

a_i = 2p', {i + 4)(i - 1)0^+2 = for i -3. (3.12) 



+ (3.13) 

where a and 7 are arbitrary differentiable functions of t. 
Note 



and 



0' = + + + + « + 4,/5'(, + + + «». (3.15) 

Substituting the above two equations into (3.7), we have: 

+107/3'(2/ + + 2i{y + + 7^(1/ + 0f. (3.16) 

Write 

/(i,y) = ^6,(i)(y + /3)\ (3.17) 

Then (3.16) becomes 

gl, . 2)(. . 1) - 21.«(, . = ^ + ^ . J^l 

+2a;7(y + /3) + 1^0 ^{y + /3)2 + 27'(y + (if + 72(y + iif . (3.18) 

Thus 2 

= foo = -«', &i = -2/5", 63 = ^, (3.19) 

64 = /?'7, ^5 = 1^, 68 = ^, (3.20) 

(i + 3)i6i+2 = for i 7^ -4,-2,-1,1,2,3,6. (3.21) 

Therefore, 

+^(!/ + + /3'7(!/ + /3)' + l^fe + + l^fe + (3.22) 
where cr and p are arbitrary functions of t. 

Theorem 3.1. We have the following solution of the equation (1.1) blowing up on 
the surface y + (3{t) = 0: 

ax 2B'x / ct^ 



+^ - «' - 2/3"(y + /?) + p(z/ + bef + ^{y + 



where a, /3, 7, a and p are arbitrary functions oft, whose derivatives appeared in the above 
exist in a certain open set o/M. 



When a — ^ — a — p — 0, the above solution becomes 



Take the trivial solution ^ = of (3.5), which is the only solution polynomial in y. 
Then (3.6) and (3.7) become 

9yy = 0, fyy = '2gt + 9^- (3.25) 

Thus 

g^ait)+my- (3.26) 

Hence 



So 



fyy ^a^ + 2a' + 2(/3' + a(3)y + (3^- (3-27) 

f^j + ay + y + ^ + j^Z/ , (3-28) 

where 7 and a are arbitrary functions of t. 

Theorem 3.2. The following is a solution of the equation (1-2): 

u^{a + Py)x + ^ + ay + y^ + ^^^^-^v' + (3-29) 

where a,j3,'~f and a are arbitrary functions oft, whose derivatives appeared in the above 
exist in a certain open set o/M. Moreover, any solution polynomial in x and y of (1-2) 
must be of the above form. The above solution is smooth (analytic) if all a, (3, 7 and a are 
smooth (analytic) functions oft. 

Remcirk 3.3. In addition to the solutions in Theorems 3.1 and 3.2, the equation (1.2) 
has the following simple solution: 

u = pXy) x^ (3.30) 
where pt^{w) is the Weierstrass's elliptic function satisfying (2.40). 

4 3-D Khokhlov-Zabolotskaya Equation 

By comparing the terms of highest degree, we find that a solution polynomial in x of the 
equation (1.3) must be of the form: 

u = f{t, y, z) + g{t, y, z)x + ^(t, y, z)x'^, (4.1) 



where f(t,y,z), g(t,y,z) and ^(t,y,z) are suitably-differentiable functions to be deter- 
mined. As (3.2)-(3.7), the equation (1.3) is equivalent to: 

^,, + ^.. = 6^^ (4.2) 

gyy + gzz - Qgi = 4^*, (4.3) 

fyy + U-m-^9t + g'. (4.4) 

First we observe that 

^ ~ {y cos a{t) + z sin a{t) + p{t)y ^^'^^ 

is a solution of the equation (4.2), where a and /3 are suitable differentiable functions of 
t. With the above ^, (4.3) becomes 

6g 8(q;'(— y since + 2;cosq;) + /3') 

^ ~ {y cos a{t) + z sin a{t)+f3{t))'' " {ycosa + zsina + ' ^ ' 

In order to solve (4.6), we change variables: 

C = cosa y + sina 2; + /3, 77 = — sin a y + cos a (4-7) 

Then 

dy = cos ad(^ — sin a drj, dz = sin a 9^ + cos a dj^. (4.8) 

Thus 

ei^j^dl^ (cos adc^- sin a dr,Y + (sin a 9^ + cos a dr,f ^ d"^ + d^. (4.9) 

Note 

dtiO = a'r^ + dt{rj) = a'(/3 - C). (4.10) 
The equation (4.6) can be rewritten as: 

gcc + gm - ^C'g = -8{a'rj + f3')C\ (4.11) 

In order to solve the above equation, we assume 

g = J2a,{t,rj)C. (4.12) 

Now (4.11) becomes 

J][((z + 2)(i + 1) - 6)a,+2 + aij = -8(a'ri + (3')C\ (4.13) 

which is equivalent to 

-4a_i + a_3^^ = -8(q;'?7 + /3'), (i + 4) (i - 1)0^+2 + = for - 3 7^ i G Z. (4.14) 
Hence 

a_i = |a_3^^ + 2(Q!'r7 + /3'), (i + 4) (i - 1)0^+2 = -a^,,^ for - 3 ^ i e Z. (4.15) 



When i = — 4 and i = 1, we get a_4^^ = oi^^ = 0. Moreover, a_2 and as can be any 
functions. 
Take 

03 = cr, a_2 = p, a_i = 2(q;'?7 + /?'), (4.16) 

oi = a_i_2j = a-2-2i = for < i e Z (4-17) 

in order to avoid infinite number of negative powers of ( in (4.12), where a and p are 
arbitrary functions of t and rj differentiable in a certain domain. By (4.15), 

(-l)^^f(c^) _ (-l)n5gf(a) 
^■=l(2^ + 5)(2^) (2^ + 5)(2^ + 3)(2fc + 1)! ' ^ " 



(-l)^gf(p) {-m2k-l){2k-3)df{p) 

= nL(2.)(2.-5) = w ■ ^'-^^^ 



Therefore, 

, = 2K, + «r' + E(-i)'l7CT^ 



fe=0 



■(2A; + 5)(2A; + 3)(2A; + 1)! 



(2A;-l)(2A;-3)af(p)C ^^^^ 



]C^'= (4.20) 



3(2A;)! 

is a solution of (4.11). 

By (4.9), (4.4) is equivalent to 

fee + Ur, - 2C-V = 2gt + g\ (4.21) 

Note 

00 

= 2{a"v + P" + {a'yPK-' - 2{a'f - 2{a'v + P'fC' + E(-l)'^''{ 

k=0 

^ y{2k + 5){2k + 3){2k + 1)! 3(2A;)! 

159f(a)C' (2A;-l)(2A;-2)(2A;-3)9f (p)C-^, , , 

^ + ^ )[ (2. + 5)(2i+l)! + 3(2.)! ' ]>" (^-22) 

For convenience of solving the equation (4.21), we denote 

2gt + g'=J2k{t,7j)C (4.23) 

i=-4 

by (4.20) and (4.22). In particular, 

6_4 = p', 6-3 - 0, (4.24) 
6_2 = 2(p, + a'/3p,) + ^, (4.25) 
= 4[a"r] + (3" + {a'f(3] - 2a'p^ + ^(a'r; + /3')p,^, (4.26) 



bo = -4(a')' + liPt^v + <^'PPmv) + + ^9%,- (4-27) 

Suppose that 

/ = 5]Q(i,77)C (4.29) 

is a solution (4.21). Then 

oo 

+ 2)(i + 1) - 2)q+2 + qJC - E ^rC, (4.30) 

ieZ r=-4 

equivalently 

(i + 3)iCi+2 = 6i - Q^^, (r + 3)rc^+2 = -Cr,,,,, r < -4 < i (4.31) 

By the above second equation , we take 

= for r < -4 (4.32) 

to avoid infinite number of negative powers of C, in (4.29). Letting i — —3, 0, we get 

6-3 = C-ZTm, bo = cor,n- (4.33) 

The first equation is naturally satisfied because c_3 = — c_5^^/10 = 0. Taking i ~ 
—2, —4 and r = — 6 in (4.31), we obtain 

Co = ^C-2r,r, " ^&-2, C_2 = ^&-4. (4.34) 



So 

Thus we get a constraint: 
equivalently, 



CO = -d',{b.,) - -6_2. (4.35) 



bo = Ja,'(&-4) - \d'^{b-2), (4.36) 



-4(q;') + -(pt^^ + a'Ppr,r,r,) + T7:9^{p)P + 



12 36 

= ^aJ(p^)-A,.-a'/3p,,,-5^. (4.37) 

Thus 

96(a,, + a'/3p,^,) + 69j(p)p + 2^^^ - 99j(p2) + I2d^^(p,r,p) = 288(a')'. (4.38) 

It can be proved by considering the terms of highest degree that any solution of (4.38) 
polynomial in rj must be of the form 

p^7o{t)+Mt)v + ^2{tW. (4.39) 

Then (4.38) becomes 

672 - 572 = 9(«')'- (4-40) 



So 



«' = |\/672 - ^ " = I 



yjei.-b'yidt, (4.41) 

where e = ±1. Replace by — /? if necessary, we can take e — 1. Under the assumption 
(4.39), 



g = pC' + 2{a'r, + /^OC"' + | + ^^(-1) 



ifc=0 



(2A; + 5)(2A; + 3)(2A; + 1)! 



and 



b-2 = 2{pt + a'Ppn) + -72P, 



Denote 



6_i = A[a"r) + (5" + (a')'/?] - 2a' p^ + -(a'77 + /3')72, 
feo = -4(a')^ + ^7^ + |. 



For any real function F{t, rf) analytic at 77 = 770, we define 



r=0 



Note 



15af (a)C='+'^ 



k=G 

= 15 



(2A; + 5)(2A; + 3)(2A; + 1)! 



{2k + h){2k)\ 



IK' 
dri dT2 



•^0 \fe=0 



{2k + 5){2k + 3){2k)\ 



(4.42) 

(4.43) 
(4.44) 
(4.45) 

(4.46) 
(4.47) 
dTi 



Jo Jo Vfe=o 



15 
T 



/ T3 / Ta / [(T(t, r] + V-iTi) + a{t, r] - ^/-lTi)]dTi dT2 dr^, 
Jo Jo Jo 



(4.48) 



159f(a)C 



fe=0 



(2A; + 5)(2A; + 3)(2A; + 1)! 



E(-i)' 



k=0 



^ i5af(a,)c^+^^ , f., ..^s^^^VK^;^^ 

(2A; + 5)(2A; + 3)(2A; + 1)! (2A; + 5)(2A; + 3)(2A; + 1)! 



^ (2A; + 5)(2A; + 3)(2A;+1)!^ '^^^^ ^ 



fc=0 



fc=0 



(2A; + 5)(2A; + 1)! 



= — C ' j j ^2 [at{t,r] + \^Ti)+at{t,r)-\^Ti)]dridT2dT3 
^ Jo Jo Jo 



2e 



} ^2 I ri[a{t, T) + a/^ti) - a{t, r] - \^Ti)]dTi dT2 
Jo Jo 



+ ^{a'ri + /?')C"' / \a{t,ri + v^n) + a{t,r]- ^/^n)]dn dT2. (4.49) 

^ Jo Jo 



Hence 



^ 1 I ^2 / [a{t,'q + ^/^Ti) + a(t,r]-^/^Ti)\dTidT2dTz, (4.50) 
Jo Jo Jo 

by (4.42) and (4.48). According to (4.23) and (4.46), we have 



*(/3,P,<x>(t,r/,C) 

225 . , / n 



K rT3 i"r2 

+ 15C-' / Ta / r2 / [(7t(i, ?7 + V^Ti) + (7t(i, 77 - V^Ti)]dTi dT2 dTs 

Jo Jo Jo 

15a'(C — 3) f'' 

+ % / / ri[a(^, ri + V^n) - a{t, rj - V^ri)]dri dT2 

s Jo Jo 

+15{a'r] + (3%-^ / r| / [a{t, rj + x/^n) + a{t, rj - V^n)]dn dT2 
Jo Jo 

(K rr?, r-T2 \ 

j '^^ j ^2 y [o{t, 7] + V^Ti) + Cr(t, 7] - \f^Ti)\dTi dT2 dTsj 

xC'{pC' + {a'v + P')C' + j)- (4.51) 

Now 2 

C-2 = ^ (4.52) 

by (4.24) and (4.34). According to (4.31) with i — —3, 0, c_i and C2 can be arbitrary. For 
convenience, we redenote 

C-i^K{t,r}), C2 = uj{t,r}). (4.53) 
Moreover, (4.24), (4.35) and (4.43) imply 

co^^-pt-a'Pp, + -^. (4.54) 
Furthermore, (4.31) and (4.44) yield 

ci = ^ - 2{a"rj + (5" + {a'f(5) + a'p, - \{a'ii + /3')72- (4.55) 
In addition, (4.31) and (4.53) gave 

-ir^df-\n) A (-l)^-(. + l)(2.)! 



2{k + 2){2k + 2)\ {k + 2){2k + 2)\ 



-ir^Sd^Hu^) ^ (-l)^-(2. + 3)(2. + l)! 



(2A; + 5)(2/c + 3)! (2A; + 5)(2fc + 3)! 



for < A; e Z. 



Set 



i=3 



oo k 



,fc=0 



39f(a;)C 



2k 



k=0 



{2k + 3)(2k + l)\ 



2k+3 



oo fe 



fc=0 i=0 



(-l)'=-^(2i + 3)(2i + 1)! 
(2A; + 5)(2fc + 3)! 



(4.58) 



Note 



^ ^ (2A; + 3)(2A; + 1)! 



fc=0 



3 
2' 



C"^ / / [a;(i,77 + v^Ti)+u;(i,?7- v^Ti)]dTidT2. 
Jo Jo 



(4.59) 



Moreover, 

*(/3,p,a)(t,r;,0) =0, 6i 



for < i e Z. 



(4.60) 



Thus 



3 K /•T'2 

^(/3,p,<T,K,a;)(*,^,C) = 2^ y ^2 y [u;(t,r? + \^ri)+a;(i,77- v^Ti)]dTidT2 



--ca^C"'[/^(^, ^ + v^Ti) + «(i, 7^ - V^n) 



oo fc 



^ ^ ^ (-i)^-^(i + i)ay-^^ar^(^(Ap,.))(t,^,o) ^,,^3 



fc=0 i=Q 

oo A; 



fc=0 i=0 



{2i + l){k + 2){2k + 2y. 

(2« + 2)(2A; + 5)(2/c + 3)! 



^ ^ ^ V -l)^--(2. + 3)gf-^af+^(vl/(,.,..))(t,r^,0) ^,,^^ 



(3.61) 



in which the summations are finite if a{t, rj) is polynomial in rj. According to (4.52)-(4.58) 
and (4.61), 

2 2 



[2{a"7j + (3" + {a'Y(3) - a' p, + -(a'r; + /3')72]C. 



(4.62) 



Theorem 4.1. In terms of the notions in (4-7), we have the following solution of the 



equation (1.3) blowing up on the hypersurface cosa{t) y + sina(t) z + P{t) = (( = 0): 

u = xV^ + [pr^ + 2(a'r^ + /5')r^ + | + ^r^ 

>^ ^3 T2 / [a{t,rj + \^Ti)+a{t,ri-V^Ti)]dTidT2dT3\x 
Jo Jo Jo 

2 2 

-[2{a"r] + P" + - a'p, + ^{a'r] + /3')72]C, (4-63) 

where the involved parametric functions p is given in (4-39), a is given in (4-41) o.i^'d P 
is any function oft. Moreover, a, k, cu are real functions in real variable t andrj, and 
^{p,p,a,K,u,){t,r]X)is given in (4.6I) via (4.51). 

When cr = = a; = 0, the above solution becomes: 

-a'(3p, + ^ - [2{a"rj + p" + (a'fp) - a'p, + '^{a'rj + /3')72]C, (4-64) 

Next we consider ^ = 0, which is the only solution polynomial in y and z of (4.2). In 
this case, (4.3) and (4.4) becomes: 

9yy + 9zz = 0: fyy + fzz = '^9t + 9'^- (4-65) 

The above first equation is classical two-dimensional Laplace equation, whose solutions 
are called harmonic functions. In order to find simpler expressions of the solutions of the 
above equations, we introduce a new notion. A complex function 

G{p) is called bar-homomorphic if G{p) — G{jl). (4.66) 

For instance, trigonometric functions, polynomials with real coefficients and elliptic func- 
tions with bar-invariant periods are bar-homomorphic functions. The extended function 
F{t, p) in (4.47) is bar-homomorphic in p. 

As (4.20), it can be proved by power series that the general solution of the ffist equation 
in (4.65) is: 

g^{a + V^p) {t, y + V^z) + {a - V^p) {t, y - V^z) , (4.67) 

where a(t, p) and p(t, p) are complex functions in real variable t and bar-homomorphic 
in complex variable p. Set 

w — y + \/^-lz, w = y — ^J^-Vz. (4.68) 



Then the Laplace operator 

dl + dl^^d^d^. (4.69) 



The second equation in (4.65) is equivalent to: 

dwdwif) ^ f + ^ = + "/-^Ptjit, w) + {ut - \^pt)it,w) 
+ ^[{a + V^p){t,w) + {a- V^p){t,w)]^. (4.70) 
Hence the general solution of the second equation in (4.65) is: 

J Wl J Wl ^ 

+ ^[((T + v^p)(t,/Ui) + (o- - ^[^p){t,-p{)f}dpx 

+(K + ^/^LJ){t,w) + {K-^/^u}){t,w), (4.71) 

where K{t,p) and uj{t,p) are complex functions in real variable t and bar-homomorphic 
in complex variable p, and Wi is a complex constant. 

Theorem 4.3. In terms of the notions in (4-67), the following is a solution polynomial 
in X of the equation (1-3): 

I'm pw ^ 

u = [{a + V^p){t,w) + {a-^/^p){t,w)]x+ _ {-[{at + V^pt)it, p,) 

+{at - V^pt){t,p:i)] + ^[((T + \^p){t,pi) + (cr - V^p){t,p:{)]'^}dpi dpi 

+ {K + \/^uj){t,w) + {K-\r-\-Uj){t,w), (4.72) 

where a{t, p), p{t,p), K{t,p) anduj{t,p) are complex functions in real variable t and bar- 
homomorphic in complex variable p (cf. (4-66). Moreover, the above solution is smooth 
(analytic) if all a, p, k anduj are smooth (analytic) functions. In particular, any solution 
of the equation (1-3) polynomial in x, y, z must be of the form (4-72) in which a, p, k and 
cu are polynomial in p. 

Remark 4.4. In addition to the solutions in Theorems 4.1 and 4.2, the equation (1.3) 
has the following simple solution: 

p,{ay + bz) x^, (4.73) 

where p,,{w) is the Weierstrass's elliptic function and a,b are real constants such that 
+ 62 = 1. 
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